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O : ABSTRACT 

O 



We estimate the expected distribution of displacements between the dark mat- 
ter and gas cores in simulated clusters. We use the MareNostrum Universe, one 
of the largest non radiative, SPH ACDM cosmological simulations. We find that 
projected 2-D displacements between dark matter and gas, equal or larger than 
Qh' the observed in the Bullet Cluster, are expected in 1% to 2% of the clusters with 

P ■ masses larger than 10^^ h~^y\.Q. The 2-D displacement distribution is roughly 

the same between redshifts < 2; < 0.5 when multiplied by a factor of (l-|-z)"^'^. 
^, ' We conclude that the separations between dark matter and gas as observed in 

the bullet cluster can be easily found in a ACDM universe. Furthermore we find 
^ ■ that the displacement distribution is not very sensitive to the normalization of 

P^ ■ the power spectrum. Upcoming surveys could extend the measurements of these 

0> ■ displacements between dark matter and gas into large samples of hundreds of 

m 



clusters, providing a potential test for ACDM. 



1. Introduction 

The standard model of structure formation describes the Universe as a three-component 
fluid (dark matter (DM), gas and radiation) plus a cosmological constant. There is a fun- 
damental intrinsic difference between the two matter components: while the DM can be 
modelled in a very good approximation as collisionless fluid, the gaseous component has 
an evolution dominated by collisional physics. As the interaction between DM an baryons 
is only through gravitational forces, one can expect a decoupling between the components 
when the collisional nature of baryons becomes important. 

From the observational point of view, the decoupling is evid ent in collisions of high mass 



clusters, observed, for the first time, in the cluster 1E0657-56 fIClowe et al.ll2006l ). dubbed 
the Bullet Cluster. For this object, a combined analysis of strong and weak gravitational 
lensing and X ray observations gives a displacement of around 200 kpc between the density 
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peak of the gaseous and dark matte r coinponents in the clus t er. More examples a r e given in 



Dupke et al.l (120071): Mahdavi et al.l (120071): I Jee et all (120071): iBradac et all (120081): IJee et al. 



( l2005al Jbl) but also see iHeymans et al.l ( l2008l ) and the peculiar case in iQin et al.l ( l2008l ). 
The standard explanation for this decoupling between DM and gas, is the collision between 
clusters of galaxies. For example, in the case of a collision between two clusters of comparable 
mass, the DM halos and galaxies go through each other without falling apart from two body 
trajectories while the intra-cluster medium, dominated by pressure effects, is stripped out of 
the DM halos. 

In the Bullet Cluster, the expected configuration is that of a massive sub-structure 
passing close to the center of the cluster, stripping away the gas from the center of both 
structures, leaving two almost gas-naked cores of dark matter and one single gas core in the 
middle. The velocity of the expected collision could be inferred obs ervationally from the 
temp erature of the bow shock t hat is usually prese n t in th e collisions ([Dopita &: Sutherlaiid 
2003h . N-body simulations (e.g JSpringel &: Farrarl J2007h . iMastropietro fc BurkertI J2008[ )) 
support the model but show also that the velocity estimation should be taken with care 
since the velocity of the encounter between dark matter structures has large statistical and 
systematic uncertainties when it is derived from the available observations. 

The Bullet Cluster has been recognized as potential stringent test for models of structure 
formation. On one hand it vindicated strongly the need for a dark matter component, but on 
the other hand it raise d doubts on the like li hood of having such extreme events in a ACDM 
universe (lNusserll2008l ). iHayashi fc White! (120061 ) attempted to calculate the probability of 
finding large velocity collisions such as the one guessed for the Bullet Cluster. They used 
a large pure DM N-body simulation (the Millennium Run, ISpringel et al.l (120051 )) and had 
to rely then on the information of the possible collisional velocities of substructures inside 
a cluster sized DM halos. Taking into account the difficulties in the interpretation of the 
data, the observed velocities can be explained within the framework of ACDM. Nevertheless, 
the frequency of such events changes from 1 to 5 in every 500 clusters, depending on the 
assumptions made on the velocities of substructures. The lowest value corresponds to the 
probability of having the sub-structure moving out of the center of the cluster. 



The downside of iHayashi fc Whitd (120061 ) is that: 



1. The mass assumed for the Bullet Cluster is lower than the most recent estimates. 
Using th e updated mass and t he same velocity constraints, the Bullet probability falls 
to 10"^ JParrar &: Rosenll2007h 



2. Their results had to be extrapolated from the data in the Millennium simulation, given 
that the abundance of clusters of comparable mass is not high enough to perform a 
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more straightforward statistical study. 



Both problems have been recently addressed by iLee fc Komatsul (|2010[ ) . They use the 
MICE dark matter only simulation with ~ 200 times larger volume than the Millennium 
Run. They use the same approach of finding large relative velocity of DM configurations. 
The large volume in the MICE simulation allows them to gather a sample of clusters with 
masses similar to the Bullet Cluster mass. Using the more stringent conditions on the 
collision velocity, they find that the probability of finding a system similar to the Bullet 
Cluster is between 3.3 x 10"^^ and 3.6 x 10"^ 

Yet, there is another indisputable and clear e xample of a Bullet- like configuration. The 
observations of the cluster MACS J0025. 4-1222 fJBradac et al.ll2008[ ) confirm the existence 
of two gas naked matter cores and one single hot gas core. As opposed to the case of the 
original Bullet Cluster, one cannot observe the X-ray emission from the bow-shock. The 
derivation of a characteristic collision velocity, from a simulation and available observations, 
can thus lead to confusing and equivocal results. In these two cases the only common and 
solid observational physical characteristic that can be measured is the separation between 
the dominant gas clump and the principal dark matter clump. 

It is expected that not all the cases of Bullet like configurati ons, will show a b ow-shock 
X-ray emission, which can constrain the sub-structure velocity f lShan et al.ll2010l ). There- 
fore, the results of the analysis of dark matter only simulations cannot properly address 
the likelihood of events similar to the cluster MACS J0025. 4-1222. If one wants to use the 
only robust and accessible observational feature, namely the distance between the domi- 
nant DM and gas clumps, to test the ACDM model predictions, then one has to resort to 
hydrodynamical simulations which include a description of the gas in clusters. 

Furthermore, as the number of observed clusters increases and will increase even faster 
with upcoming surveys like eRosita, LSST and Pan-STARRS, we are approaching an era in 
which a statistical study of this decoupling would become possible through weak lensing and 
galaxy number counts as a proxy for DM distribut ion and the X-ray emission as a tracer of 
the gas component. As a first step in this direction IShan et al.l (120 lOf ). studied a compilation 
of 38 clusters. They found 13 objects in which the separation between the center of mass 
of the gas clump and of the dark matter halo is greater than 50 kpc and 3 more clusters in 
which the separation is greater than 200 kpc. 

The purpose of this paper is to aboard the likelihood of this decoupling from the point of 
view of the distance between the DM and gas cores in clusters measured in one of the largest 
volume hydrodynarnical s imulations done up to now: the so-called MareNostrum Universe 
dCottlober fc Yepesll2007l ). 
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The structure of the paper is as follows: In Section [2] we present the main features of 
the simulations, in Section |3] we define the algorithm to find the separation between gas 
and DM structures. In Section H] we present the probability of finding a cluster with a 
given projected displacement on the sky. We inspect in Section [S] the dependence of that 
probability with respect to the normalization of the initial power spectrum. In Section [6] we 
scrutinize whether the Bullet Cluster should be considered as a challenge to A CDM in the 
context of our study. We summarize our results and state our conclusions in Section [3 



2. Simulation 



The analysis presented in thi s paper uses mainly the M arenostrum Universe SPH cos- 
mo logical simulation described in iGottlober fc Yeped (120071 ) . Thi s non-radiative simulation 
was run using the GADGET2 TreePM-|-SPH code ISpringell (120051 ) and follows the evolution 
of gas and dark matter from z=40 to z=0 in a comoving cube of 500 h~^Mpc on a side. The 
cosmology used corresponds to the spatially flat concordance model with the following pa- 
rameters: the total mass density Qm = 0.3, the baryon density Qb = 0.045, the cosmological 
constant Q^ = 0.7, the Hubble parameter h = 0.7, the slope of the initial power spectrum 
n = 1 and the normalization erg = 0.9. The number of particles used for each of the DM and 
gas component was 1024^, resulting in a mass of 8.3 x 10^/i~^Mq for the DM particles and 
1.5 X 10^/i~^Mq for the gas particles. The spatial force resolution was set to an equivalent 
Plummer gravitational softening of 15 /i~^kpc comoving. 

In order to study resolution effects on the distance determination, we have also used 
a lower resolution simulation run from the same initial conditions as the MareNostrum 
Universe. Here, the number of particles is 2 x 512^ and the gravitational smoothing length 
is correspondingly e ~ 25 h~^kpc . We have performed two more simulations with the same 
number of particles but different cosr nological parameter s (equivalent to those obtained from 
fitting the 3rd year of WMAP data (jSpergel et al.l 120071 )) and two different power spectrum 
normalizations, a low erg = 0.75 and a medium erg = 0.8. The objective is to quantify 
the effect of that cosmological parameter on the separation probabilities. In Table [1] we 
summarize the main features of our simulations. More details of these simulations can be 



found elsewhere ( JYepes et al.l 120071 ). 



We analyse in the following a non-radiative simulation of a large cosmological volume. 
It is well known that this approximation fails if one considers the inner regions of clusters 
(within about 0.1 of the virial radius). However, the outer regions of clusters are nearly 
self-similar and different simulations agree in a remarkably successful description of the tem- 
perature profile in the outside regions of clusters. There are many studies of the formation 
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of isolated clusters which take into account cooling processes and star formation, however 
these simulations still fail to describe the inner structure of clusters correctly. Most probably 
additional processes, hke AGN feedbacks, should be taken into account for a proper descrip- 
tion. In a recent review, Borgani and Kravtsov (JBorgani fc Kravtsovll2009l ) discuss in great 
detail the issue of cosmological simulations of galaxy clusters. In any case, the statistical 
study of cluster mergers needs a large simulation volume and large number of particles to 
properly resolve the merger physics. This task is not yet possible to be done when radiative 
physics is also considered. Furthermore, there are studies that suggest that the incomplete 
and rough non-radiative des cription provides a fair approximation to the dynamics of the 
gas peaks (jPoole et al.ll2006l ). 



3. Separation Measurement Algorithm 

The aim of our analysis is to obtain the distribution of displacements between the domi- 
nant density peaks of the DM and gas components of each selected cluster in the simulation. 
We describe first the algorithm used on a set of hierarchical friend-of-friends (FOF) catalogs. 
Then, we describe briefly a different and more complex algorithm of structure detection, 
which, nevertheless, gives the same statistical results than the hierarchical FOF. 



3.1. FOF Based Algorithm 



We have analys ed all simulations with the hierarchical friend-of-friends (FOF) algorithm 
( iKlypin et al.lll999l ) used for the DM and gas particles, separately. The basic FOF groups 
where identified using a linking length of 6 = 0.17 times the mean inter particle separation 
both for the gas and the DM distribution. The substructures inside the FOF groups are 
defined as FOF objects constructed with shorter linking lengths 6.„ = 6/2" (n=l,2,3). We 
name these catalogs FOF^^^ and FOF^"* with (n=0,l,2,3) for the DM and gas components 



Name 


Particles 


0"8 


VIm 


^A 


n 


MN-1024 


2 X 10243 


0.9 


0.30 


0.70 


1 


MN-512 


2 X 512=^ 


0.9 


0.30 


0.70 


1 


MN-512-MS 


2 X 512^ 


0.8 


0.24 


0.73 


0.95 


MN-512-LS 


2 X 512^ 


0.75 


0.24 


0.73 


0.95 



Table 1: Summary of the main characteristics of the four simulations used in this work. 
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Redshift # Selected Clusters Most massive cluster 
z (IQi^ /i-^Mq) 

25 
20 
13 

Table 2: Analyzed redshifts z, number of selected clusters with M > 10^^ h~^MQaiid mass 
of the most massive cluster at the given redshift in the simulation volume 






4063 


0.3 


2662 


0.5 


1826 



Mass bin Minimum Mass Maximum Mass # Halos # Halos 7^ Halos 

Mh[l0^^h-^MQ] Mh[10^^h-^MQ] z = z = 0.3 z = 0.5 



1 
2 



1.00 
2.00 



1.12 
25.0 



616 
1403 



463 
759 



352 
409 



Table 3: Mass bins ranges for the different redshift analysis. The last three columns show 
the number of halos in that bin at a given redshift. 





Mass- 1.0 lO^^h-iMsun z = 0.3 




360.8 








:/y^;^!C^^^^^^^>V\gas f. 


71.8 


^60.6 

Q. 

5 




1 .T ~/-v "■ \ 


|71.6 

T 

N 

71.4 


f360.4 




■ i~\- ^ r 


360.2 




71.2 


63.0 63.2 63.4 63.6 






x[h-iMpc 







Mass- LOIOi^h-ilVlsun z = 0.3 




li I I I . I ■'. t:rC I . ^r . lZ, I ; 



63.0 63.2 63.4 63.6 

X [h^^Mpc] 



Fig. 1. — Linearly spaced surface density contours for the DM (dashed blue) and gas (con- 
tinuous red) components of one of the halos with a high displacement seen in two different 
projections. The arrows are proportional to the velocities respect to the center of the gas 
core. The cluster of mass 1.0 x 10^"^ /i^^MqIs located at redshift z = 0.3, the coordinates 
are comoving and the relative speed between the two dark matter cores is 560 km/s. In this 
case, the physical displacement between the centers of gas and DM is of the order of the 
observed cluster MACS J0025.4-1222. 
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respectively. For more details on the simulation, the FOF identification and general cluster 



properties we refer the reader to iGottlober fc YepesI (120071 ) 



For our analysis we have selected all halos with masses larger than 10^^h~^MQ, at four 
different redshifts, z= 0, 0.3 and 0.5. The number of selected halos for every redshift, and 
the highest halo mass are quoted in Table [2l 

Our objective is to measure the physical separation between the dominant gas clump in 
the cluster, with respect to its predominant DM structure. We find the separation based on 
the FOF catalogs described above. The main advantage of the FOF algorithm is that it can 
be easily applied to the subsets of DM particles and gas particles and identifies gas clumps 
even if they don't have a pronounced density peak. Moreover, the set of linking lengths 
in the hierarchical FOF allows to find objects at different overdensities. In order to test 
the robustness of the results from the algorithm based on FOF catalogs we have compared 
them with those obtained using a spherical overdensity algorithm, the Amiga Halo Finder, 
described in the next section. 

We define the clusters as the objects detected in the DM FOFq catalog. FOFq provides 
objects with about virial overdensity, whereas FOFj provides substructures with about 8* 
higher densities. In the MN-1024 simulation a typical object of lO^^/i~^M0 at the low mass 
end of our sample consists of ~ 10000 DM particles and about the same number of gas 
particles. In the series of low resolution simulations we consider only clusters with masses 
larger than 3.2 x 10^^/i~^Mq which consists of ~ 4000 DM particles and about the same 
number of gas particles. 

In all cases we find that the clusters exhibit a massive dominant gas clump and various 
DM clumps. In both DM and gas catalogs we define the dominant clump position in the 
following way. We identify first the center and radius of the cluster as given by the DM 
FOFq catalog. Next we find the most massive group in the FOFi catalog with its center 
inside the radius defined previously. The clump defines a new center and radius. We iterate 
this procedure until we cannot find a new group in the following FOFj catalog. The outputs 
of the algorithm are the mass and radius of the cluster defined from the DM FOFq catalog 
and the centers of mass of the dominant DM and gas clumps as well as their separation. 

It is well known that the FOF algorithm at the level FOFq will connect two clusters 
that are starting to merge as one FOF object. Since such cases could lead to spurious large 
separations between the DM and gas centers they have been excluded from the analysis. 
These cases represent always less than 3% of the total number of clusters. 



3.2. AHF Based Algorithm 

As an example of a spherical overdensity halo finder we have used the Amiga Halo 
Finder (AHF) u which identifies both halos and subhaloes. AHF locates local overdensities 
in a adaptively smoothed grid as prospective halo centers. AHF is described in detail in 



KnoUmann fc Knebd (120091 ) 



The AHF algorithm proceeds by creating a hierarchy of grids. High refined grids are 
spanned in high density regions until a minimum number of particles per grid-cell is reached. 
The algorithm includes a criterion to span a new mesh refinement, which is usually expressed 
as a density threshold. To find the positions of the gas and DM clumps we run AHF over 
the gas and dark matter distribution separately. Thus the positions can be associated to the 
regions where the corresponding mesh refinement reached its highest level. 

In the studied range of separations the two algorithms give the same answers within 
the simulation uncertainties. We are confident that our results are robust for the range of 
observed Bullet Cluster displacements. We proceed the rest of the analysis in this paper 
using the hierarchical FOF algorithm, which is less computationally expensive and much 
faster than the AHF one. 



4. Cumulative Probability Distributions 

We found simulated clusters that show larger separations than the observed Bullet Clus- 
ter or the cluster MACS J0025. 4-1222. We checked by visual inspection all the cases where 
the separation between DM and gas peaks is larger than 200 kpc. All of them correspond 
to perturbed clusters. In Figure [1] we show one example. The cluster was selected because 
of its striking similarity to the observed configuration in MACS J0025. 4-1222. 

We have also checked the dynamical state of two different samples of the clusters at 
redshift zero. The first sample included those with large displacements (> 200 kpc) and the 
second sample comprises all clusters with low displacement (< 10 kpc) between gas and dark 
matter. We found that, indeed all clusters with large displacements come from plunging 
substructures. Moreover, those clusters which show the most extreme displacements also 
present a rich variety of merging configurations, with triple DM cores in some cases. On the 
other hand, in all the cases of clusters with low displaments, the dominant gas blob has a 



^It is a MPI+OpenMP hybrid halo finder to be downloaded freely from 



http: //www.popia.ft .uam. es/AMIGA 
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unique dark matter core associated to it, with practically zero relative velocity between them. 
A detailed analysis and classification of all configurations giving rise to large displacements 
in the simulated cluster will be presented in a forthcoming paper. 

From the simulations, we know the three dimensional positions of the dominant DM 
and gas clumps. In order to compare our results with observations, we use 2D projections of 
the 3D displacements into 10 different randomly selected directions. We quantify the results 
using the cumulative distribution of the projected 2D displacements. Using the 10 different 
2D displacement distributions we bin logarithmically the data in distance and find the mean 
value and dispersion for the corresponding cumulative distribution. 

In the following subsections we study the scaling with mass and redshift of the calculated 
2-D displacement cumulative distributions. 



4.1. Scaling with Cluster Mass 

In Figure |2] we present the cumulative distribution of projected 2D displacements in the 
MN-1024 simulation. We bin the selected clusters in mass, and calculate the distribution 
for the 2D physical displacements. The two different mass bins are described in Table (E])- 
The mass bins are selected in such a way as to have a similar number of clusters in each bin 
at z = 0.5. We find that, at ^ = 0, the fraction of clusters with small physical separations 
(< 200kpc) increases with cluster mass. In constrast, at redshifts z = 0.3 and z = 0.5 we 
cannot assert that this trend with halo mass exists at all. In any case, in the range of the 
Bullets displacements, the distributions of different mass bins are less than 1% apart for all 
the redshifts considered here. As can be seen in the figure, the fraction of simulated clusters 
with a displacement equal or larger than the observed ones is always between 1% and 2%. 
The lowest value corresponds to the less massive clusters at z = and the highest value is 
associated with clusters at high redshifts z ~ 0.5. 

The statistical properties of the displacements at redshift z = 0.3, for masses larger 
to 1.0 X lO^^/i~^M0, is nevertheless very poor. We find a small sample of 8 clusters. An 
insufficient number to make a solid statistical statement. In spite of that, we found 1 cluster, 
from the sample of 8 clusters, with a 3D displacement larger than the one measured in the 
Bullet Cluster. 
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Fig. 2. — Cumulative distribution of 2D displacements between dark matter and gas peaks 
at redshifts z = 0, 0.3 and 0.5 as indicated in the corner of each plot. Triangles represent 
cluster more massive than 2.00 x 10^^h~^MQ, squares are clusters with masses between 
1.00 X lO^^h-^MQ and 1.12 x IO^^/i-^Mq. Mh is given in units of 1 x 10^^ /i-^Mq. The 
upper row shows the 2D displacements in physical units, the lower row normalized by the 
virial radius of the cluster. The error bars associated to each point are calculated from the 
variances of 10 different 2D cumulative distributions, computed for different projections of 
the 3D separations measured in the simulation. The vertical lines in the upper row indicate 
the displacements for the observed Bullet Cluster (Bullet- A) and the cluster MACS J0025.4- 
1222 (Bullet-B). The occurrence of a similar or larger displacement is then to be expected 
in between 1-2 clusters out of 100. 
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Fig. 3. — Cumulative distribution of 2D physical displacements between dark matter and 
gas peaks as a function of redshift scaled by \/\ + z: z = (hexagons), 0.3 (triangles) and 
0.5 (squares). All halos more massive than lO^^/i~^M0 were taken into account to make this 
plot. The vertical lines indicate the measured displacements for the observed Bullet Cluster 
(Bullet-A) and the cluster MACS J0025.4-1222 (Bullet-B). The red dashed line corresponds 
to the fitting formula given in Eq.([T]). 
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4.2. Scaling with Redshift 

In Figure [3] we show the redshift evolution of the cumulative distribution of 2D physical 
displacements. All clusters more massive than 10^^/i~^MQwere taken into account for the 
construction of these distributions. As can be seen in this Figure, when we consider the 
rescaled variable d/^Jl + z, where z corresponds to the redshift, the cumulative distributions 
turn out to be the same for all the redshifts analyzed. 

The shape of the cumulative distribution, seems to follow a power law of the displace- 
ment until it reaches a critical value, close to the Bullets separations ~ 300 kpc, and then 
falls off exponentially. A Schechter-like function turns out to be a good fit to the numerical 
results. In terms of the redshift scaled variable dz = d/\^l + z: 

'dzY f d 
d- ^"P "Z 



P, 



2D 



Pn 



1) 



with parameters Pq = 0.04, d^, = 200 kpc and a 
Figure [31 in the range 100 < dz/kpc < 500. 



-1.0, fit very well the data shown in 



4.3. A toy model to explain the separation dependences 

The redshift scaling of the cumulative probability of displacements shown in Fig. [3] calls 
for an explanation. To this end we resort to previous numerical studies that have followed the 
dynamics of colliding clusters, in a general context rather than just focusing on reproducing 
the Bullet Cluster. 



The are two important studies that we take as basic references to this end. iTormen et al. 



( 120041 ) follow ed the prop e rties o f cluster satellites from hydrodynamical cosmological sim- 
ulations and iPoole et al.l ( 120061 ) studied the impact of mergers on relaxed X-ray clusters, 
focusing on the dynamical evolution and the emergent transient structures. Both studies 
were based on numerical simulations which include gas dyna mics, but only the l atter in- 
cluded cooling and star formation as well. While the work of iTormen et al.l (l2004r) is done 
in an explicit cosmological setup, in the numerical simulations of IPoole et al.l ( 20061 ) the cos- 
mological context is taken into account in an implicit m anner through an a dequate selection 
of the initial conditions. It is relevant to highlight that iPoole et al.l (120061 ) fi nd that the re 
suits f rom their simulations compare well with the simple model predictions in lTormen et al. 
( I2OO4J ). despite the large differences in the modelling of the baryonic physics. 



More specifically, iPoole et al.l ( 120061 ) find t hat the disrpu t ion of the gas substructure is 
quantitative similar to the results obtained by iTormen et al.l ( l2004l ) even though the effect 
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of cooling allows for the formation of denser resilient gas cores. We should also expect 
that qualitatively, the results presented in this paper would not change much if cooling was 
included in our simulation. However, a quantitative statement on how the displacements we 
measured are affected by cooling and star formation processes, would require to run new, 
high-resolution re-simulations of our clusters. This analysis is well beyond the scope of the 
present work and will be devoted in further publications on the same topic. 



In most of the cases studied by lPoole et al.l (120061 ) the central ICM of the merger remnant 
oscillates with respect to the dark matter distribution. The oscillation of the main ICM 
remnant can be seen exclusively between the first core interaction and the relaxation of the 
cluster, a process which lasts between IGyr to 3Gyr depending on the mass ratio and the 
impact parameter of the collision. During this time interval, the core oscillation presents 
only one large amplitude separation, with the subsequent separation being around half of 
the first one. A third oscillation shows always an even shorter separation on the order of 50 
kpc. One could then simplify the dynamics of the gas-DM separation with the description 
of damped oscillations, characterized by an initial amplitude, an intrinsic frequency and a 
damping timescale. Based on the above results, we try now some simple estimations of the 
scalings of the physical magnitudes in the model. 



Let us consider a momentum conserving collision between two gas blobs of mass Mg 
and rrig {Mg > rrig), with initial velocities Vg = and Vg > 0, respectively. Assuming a 
totally inelastic collision, the velocity of the final blob of mass Mg + rUg can be expressed 
a.s Vi = Vgmg/{Mg + mg). With the mass ratio definition x = Mg/rrig, we can rewrite 

Vi=Vg/{X+l). 

N-body simulations of halo mergers have shown that the average velocity of an accreting 
system into rich cluster type halos is of the same order than the circular velocity, K, of th e 
main cluster, independent of its mass (JTormerull997l : IVitvitska et al.ll2002l : iPoole et al.ll2006l ). 
We can then approximate Vi ~ Vc/{x + 1). Thus, the first initial amplitude of the DM core- 
core interaction A oc Vc/{x + 1). 

The intrinsic frequency of the core-core oscillation can be estimated to be proportional 
to the apocentrico time for the dark matter component of the satellite. The apocentric time 
encodes the information on t he typical dynamical time-scale around the center of the main 
cluster. iTormen et al.l (j2004j ) find that the apocentric times depend only on the mass ratio 



X. 



^Thc apoccntcr is identified as the first relative maximum in orbital distance measured after the first 
relative minimum in the satellite orbital distance. 
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The last element in this description, the intrinsic damping time-scale can be physically 
related to the same process that un binds gas and dark matter from the satellite and force 



their merging. iTormen et al.l ( 120041 ) find that the decay time to unbind self-bound gas and 



dark matter particles from the satellite depends mainly on the mass ratio x. 

Summarizing, the separation between gas and DM depends primarily on the circular 
velocity Vc of the main cluster and the mass fraction x of the merging systems. Thus, the 
larger amplitudes are in principle associated to clusters with larger circular velocities, but 
depend also on the distribution of mass ratios x. This might be the reason why we do not 
see a clear dependence with cluster mass at all redshifts. 

Moreover, there is also a redshift dependent scaling relationship between the circular 
velocity Vc in physical i nagnitudes, and the halo mass Mh. The redshift scaling is such that 



docVcOC M^-^^y/lT^ flTormen et al.ll2004h . 



In Section I4l2] we considered the displacements of all halos. In this case, the cumulative 
probability distribution is dominated by the low mass halos simply because of its higher 
abundance, meaning that we are probing a similar mass range ~ 10^'^/i~^MQat different 
redshifts. From the numerical data we find that the distribution of d/\/l + z is redshift 
independent. This redshift scaling is consistent with the one derived from our toy model. 



4.4. How many Bullets should we expect? 

One can now ask the following question: which is the fraction of clusters that will show 
a separation larger than a given distance? If the volume limited sample is sensitive to halos 
down to 10^^ h~^MQ, the answer to this question is given by Figure [S] a fraction between 
1% and 2% of the clusters presents a separation equal or larger than the Bullet, depending 
mostly on the redshift, and weakly on the cluster mass, as seen in Figure [2J If we extrapolate 
these results to clusters with masses larger than in our simulations, we would expect that 
the fraction of bullet-like displacements will not be very much different (within a factor of 
~ 2) of the 1-2 % shown here. 



5. Dependence on Resolution and Power Spectrum Normalization 

We will quantify now the dependence of our results on the numerical resolution of 
the simulations. To this end we compare the cumulative distribution of 2D displacements 
of DM and gas centers at redshift 2; = of the simulations MN-1024 and MN-512. The 
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only difference between the two simulations is the number of particles and the gravitational 
smoothing. The same seed has been used for the initial conditions. 

For this comparison we use the highest mass bin in which we have, in the low resolution 
simulation, about 4000 DM particles per cluster. This is enough to determine the density 
peaks with high precision. In the left panel of Figure H] we show that there is a good 
agreement between the cumulative probability distribution of the DM-gas displacements in 
the two simulations with different resolution. 

Therefore, we are convinced that our results are not likely to be affected by numerical 
resolution as long as we have 10.000 or more gas and DM particles available in the clusters 
under consideration. 

Now, we want to use the low resolution simulations to study the dependence of our 
predictions on the underlying cosmological model. Here we focus again on clusters with more 
than 4000 DM and gas particles. We compare three models with different normalisation and 
matter co ntent which are ro ughly compatible with the WMAPl, 3, and 5 predictions (see 
Tabled]). lYepes et al.l (120071 ) have shown that different ag give rise to strong variations in 
X-ray cluster abundance. 

We repeat the analysis presented above using the low resolution simulations MN-512-LS, 
MN-512 and MN-512-MS. In the right panel of Fig. IHwe show the 2D probability separation 
distributions for the three lower resolution simulations at redshift z = 0. We cannot argue 
for a significant difference in the probability distributions derived from simulations with 
different ag . Although the total number of cluster is different in each simulation, the fraction 
of clusters showing a give displacement of DM-gas peaks remains almost unchanged. 



6. Are Bullet Clusters a Challenge to ACDM? 



The recent work by lLee fc Komatsul (120 lOl ) claims again the difficulty represented by the 
Bullet Cluster to the AC D M cosmology. Their an alysis is very similar in spirit to the work of 
Hayashi fc White! ( l2006l ). iLee fc Komatsul (120101 ) analyze a large volume dark matter only 
simulation looking for high mass clusters that were not found in the relatively small volume 
spanned by the Millennium Simulation. Th e relative velocity of the s tructu re colliding with 
the main cluster is taken from the study of iMastropietro &: BurkertI ( l2008l ). 



The study undertaken by IMastropietro fc BurkertI (120081 ) explores the parameter space 
in a set of idealized cluster collisions. The best mass ratio, impact parameter, and relative 
velocities of the collision are constrained by assuming that a set of four properties is reason- 
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ably reproduced: i) the displacements between the X-ray peaks and the mass distribution, ii) 
the morphology of the bow shock, iii) brightness of the shock and iv) the morphology of the 
main gas peak. Most of the simulations were performed with a non-radiative treatment of 
the gas physics. Two different runs out of the total 13 included a crude implementation for 
cooling. Leaving aside the well known overcooling problem, which affects the main gas peak, 
the brightness of the shock is strongly modified. Thus the simulated properties of the bow 
shock are still not well understood and are not yet robustly described by these simulations. 

One of the most robust simulated properties is related to the main gas peak. It is 
well known that a model of cooling without inclusion of energ etic feedback does not giy e 



a proper description of the gas in the inner part of a cluster f lBorgani fc Kravtsovll2009l ). 



Nevertheless, already the non-radiative d escription of the g as is a reasonable approximation 



for the dynamics of the main gas peak (jPoole et al.ll2006l ). A complete description of the 



physics of shocks, should include cooling, star formation and feedback. The fact that the 
bow shock properties are not robustly described neither in adiabatic nor in radiative runs, 
weakens the case for its selection as one of the criteria to constrain the success of a certain 
cosmological model in reproducing the observations. We doubt that the disagreement of the 
predicted and observed properties of colliding clusters depends only on the parameters of 
the collisions which are determined by the underlying ACDM model. 

In other words, simulations of structure formation in the ACDM paradigm still do not 
predict robustly the present very detailed X-ray observations of colliding clusters. Therefore, 
at present it it is difficult (or impossible) to decide whether or not the complete set of observed 
properties of the cluster 1E0657 — 56 show a significant anomaly. 

In this paper we have taken a much more modest objective, of selecting a single robust 
property in ACDM simulations of clusters, namely the distance between the main dark 
matter and gas peaks, to work out the significance of the observed anomalies in clusters like 
1E0657— 56 and MACS J0025.4-1222 . We find that large DM-gas displacements are fairly 
common in a ACDM cosmology. Furthermore our theoretical estimates of the fraction of 
cl usters with large d isplacements are roughly consistent with the recent observational study 



of IShan et al.l (120 lOf ). Therefore, we consider that Bullet Clusters should not be considered 



as a challenge to ACDM. 



7. Summary and Conclusions 

Inspired by the displacement between gaseous and dark component observed in a num- 
ber of clusters of galaxies, we calculate the distribution of these displacements in a large 
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cosmological SPH simulation including dark matter and gas. 

Theoretically, the mechanism giving rise to these displacements involves a high velocity 
substructure passing through the center of the parent cluster halo. Assessing the feasibility 
to accommodate that kind of events in the ACDM context, could be used as a test for the 
validity of the model. Unfortunately, deriving the relative velocity from observations is a 
non-trivial task, because of statistical and systematical uncertainties, making difficult the 
comparison to theoretical models. 

The separation between the DM and gas components is a much better defined quan- 
tity, mostly affected by statistical errors in the measurement process, but not so much by 
systematic or model-dependent uncertainties. 

In this paper we derive, for the first time, the expected physical 2D separation distri- 
bution for clusters with a DM mass larger than 10^^/i~^Mq. We find that around 1% to 2% 
of these clusters show DM-gas separations equal or larger than the observed in the Bullet 
Cluster. Thus, the existence of Bullet Clusters should not be considered as a challenge to 
the ACDM model as it has been recently claimed. 



Even though the fract ions of expected bu 



lets i n our ACDM simulation basically coincide 



with the previous work of iHayashi fc White! (120061 ). the comparison between the two results 



should be taken with a bit of caution. As we mentioned before, the definition of a Bullet 



cluster in iHayashi fc White! (12006! ) is different than in our case. Moreover, the dynamical 
origin of the large displacements in our clusters seems to be more varied than the scenario 
of a single dominant substructure passing through the cluster. Nevertheless, the striking 
coincidence of the predicted bullet fractions between the Millenium N-body simulation and 
the MareNostrum SPH simulation suggests that indeed it is possible to produce a bullet-like 
configuration with substructure velocities as low as ~ 2400 km/s. 

The a pproach we use and the r esult s we found are above a ll a necessary complement to 
the work of IHayashi fc White! (12006! ) and JLee fc Komatsul (|2010! ). We give information about 
the configuration space rather than the velocity space, which can be compare d in a more 
straig htforward way to observations as it has been done in the recent work of IShan et al. 
(120 lOl ) on a small sample of 38 clusters, who find an acceptable agreement with our theoretical 
predictions. The full distribution of these displacements can be considered as new prediction 
of the ACDM model, which could eventually be compared against observations. Upcoming 
large optical and X-rays surveys can make this feasible in the next years. 
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Fig. 4. — Left panel: Cumulative distributions of 2D displacements for simulations with 
different mass resolutions. Righ panel: Cumulative distributions of 2D displacements for 
simulations with different power spectrum normalizations ag- Mh is given in lO^^/i~^M0. 



